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Abstract This paper investigates the connection between discrete and continuous models 
describing prion proliferation. The scaling parameters are interpreted on biological grounds 
and we establish rigorous convergence statements. We also discuss, based on the asymptotic 
analysis, relevant boundary conditions that can be used to complete the continuous model. 
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Introduction 

| The modelling of intracellular prion infection has been dramatically improved in the past 

£Nj ■ few years according to recent progress in molecular biology of this pathology. Relevant 

models have been designed to investigate the conversion of the normal monomeric form of 
£ — • the protein (denoted PrPc) into the infectious polymeric form (denoted PrPsc) according to 

<^ ' the auto- catalytic process: 

CT\ PrPc + PrPsc — > 2PrPsc, 



O 



>< 



in fibrillar aggregation of the protein. These models are based on linear growth of PrPsc 
polymers via an autocatalytic process 



The seminal paper by Masel et al. 



Eig96 



MJN99] proposed a discrete model where the prion 



population is described by its distribution with respect to the size of polymer aggregates. 
The model is an infinite-dimensional system of Ordinary Differential Equations, taking into 
account nucleated transconformation and polymerization, fragmentation and degradation 
of the polymers, as well as production of PrPc by the cells. This model consists in an 
aggregation fragmentation discrete model. In full generality, it writes as follows: 



d °° 

_ = A - 71; - v ^ nut + 2 ^ ^2 ik i,j^u h 



j>n i<n (T\ 

dt 

j>l 

(with the convention t„ _iu„ _i = 0). Here v represents the quantity of healthy monomers 
(PrPc), Ui the quantity of infectious polymers (PrPsc) of size i, i.e. formed by the fibrillar 
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aggregation of i monomers. We thus have i ^ no > 2, where no represents the minimal size 
for polymers: smaller polymers are considered to be unstable and arc immediately degraded 
into monomers, as the last term of Equation {T]) for v expresses. 7 and /ij are the degradation 
rates respectively of monomers and polymers of size i. A is a source term: the basal synthesis 
rate of PrPc. (3i is the fragmentation rate of a polymer of size i, and kjj is the repartition 
function for a polymer of size i dividing into two polymers of smaller sizes j and i — j. Finally. 
vTj is the aggregation speed of polymers of size i, which is supposed to depend both on the 
available quantity of monomers v and on a specific aggregation ability Tj of polymers of size 
i. 

In the original model MJN99 , the degradation rate of polymers ^ and the aggregation 
rate Tj were assumed to be independent of the size i, the fragmentation rate satisfied f3j = 
/3{j — l) for a constant (3 and ki,j was a uniform repartition over {1, . . .j — 1}, i.e., kij = jrj- 
These laws express that all polymers behave in the same way, and that any joint point 
of any polymer has the same probability to break. It allowed the authors to close the 
system into an ODE system of three equations, which is quite simple to analyze. However, 
following recent experimental results such as in |SRH + 05] . and their mathematical analysis 
in [CLO + 09l ICLD + 08] , we prefer here to consider variable coefficients in their full generality. 
Following the ideas of |DC95j . we can consider, under reasonable growth assumptions on the 
coefficients, the so called admissible solutions, i.e., solutions obtained by taking the limit of 
truncated systems (see Appendix IA.3|) . 

Recent work by Greer et al. analyzed this process in a continuous setting [GPMW06] . 
They proposed a Partial Differential Equation to render out the above-mentioned polymer- 
ization/fragmentation process. It writes 



— = A - jV - V j 



r(x)U{t, x) dx + 2 / / 

xq J x—xq J y— 



dU d [°° 

— = -(,(x)U(t,x)-p(x)U(t,x)-V—(TU)+2j 0(y)k(x,y)U(t,y)dy. (3) 

The coefficients of the continuous model ©© have the same meaning than those of the 
discrete one ([T]); however, some questions about their scaling remain, and in particular about 
the exact biological interpretation of the variable x. 

The aim of this article is to investigate the link between System (TTJ) and System ©(HI). 
We discuss in details the convenient mathematical assumptions under which we can ensure 
that the continuous system is the limit of the discrete one and we establish rigorously the 
convergence statement. We also want to discuss possible biological interpretations of our 
asymptotic analysis, and see how our work can help to define a proper boundary condition 
at x — xo for System JS])©. Indeed, Eq. §3§ holds in the domain x > xq and, due to the 
convection term, at least when V(£)t(xo) > a boundary condition is necessary to complete 
the problem. 

In Section [T] we first recall general properties and previous results on the considered 
equations. In Section [21 we rescale the equations in order to make a small parameter s 
appear, and state the main result: the asymptotic convergence of the rescalcd discrete system 
towards the continuous equations. Section [3] is devoted to its proof, based on moments a 
priori estimates. Sections U and [5] discuss how these results can be interpreted on physical 
grounds. We also comment the issue of the boundary condition for the continuous model. 
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1 Basic properties of the equations 

All the considered coefficients are nonnegative. We need some structural hypothesis on k 
and kj i to make sense. Obviously, the hypothesis take into account that a polymer can only 
break into smaller pieces. We also impose symmetry since a given polymer of size y breaks 
equally into two polymers of size respectively x and y — x. Summarizing, we have 

ki,j>0, k(x,y)>0, 

kj t i = for j > i k(x, y) = for x > y, 

k jt i = ki- jt i, k(x,y) = k(y - x,y), (5) 

■J'- 1 rV 

^ kjj = 1, / k(x,y)dx = 1. (6) 
i=i J ° 

(Note that (0| and © imply that < fc^.j < 1.) Classically, these two conditions lead to a 
third one, expressing mass conservation through the fragmentation process: 

2 ^ ik hj = i> 2 xk(x,y)dx = y. (7) 
i=i J ° 

The discrete equation belongs to the family of coagulation-fragmentation models (see 
|BC90j . [BCP86] ). Adapting the work of [BC901 IBCP86] to this system, we obtain the 
following result. It is not optimal but sufficient for our study. 

Theorem 1 Let ki.j satisfy Assumptions We assume the following growth assump- 

tions on the coefficients 



There exists K > 0, a > 0, m > and < 9 < 1 such that 
< ft < Ki a < m < Ki m , 0<Ti< Ki e . 

The initial data v > 0, u® > satisfies, for a = max(l + m, 1 + 6, a) 

oo 

i a u° t < +oo. 

i=n 

Then there exists a unique global solution to {!]) which satisfies for all t > 

OO OO n £ pi OO 

v(t) + iu{(t) = v° + iu® + \t — / r yv(s) ds — / ifiiUi(s) ds. 



(8) 



(9) 



A sketch of the proof is given in Appendix IA.3I Let us introduce the quantity 

oo 

p(t) =v{t)+^iUi{t), (10) 

n 

which is the total number of monomers in the population. Equation is a mass balance 
equation, which can be written as 

^ oo 

— p = x-^v(t) - ^2 mu>i{t). (ii) 

i=n 

Similarly for the continuous model we define 

g{t) = V{t) + / xU(t,x)dx. 
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The analogue of (fTTjl would be 



g(t) - q(0) = Xt - / jV(s)ds- / xn(x)U(t,x) dx. (12) 

Jo Jo Jx 

In fact, the argument to deduce (JT2J) from the system (JH)© is two-fold: it relies both 
on the boundary condition on {x = xq} for ([3]) and on the integrability properties of the 
fragmentation term 

x X (2 f f3(y)k(x, y)U(t, y) dy - p{x)U{t, x) 

the integral of which has to be combined to © by virtue of ([7]). The question is actually 
quite deep, as it is already revealed by the case where fj, = 0, r = and xq = 0. In this 
situation it can be shown that ([3]) admits solutions that do not satisfy the conservation law: 
Jg 00 xll(t, x) dx = J °° xU(0, x) dx, see [DS96 . Hence, (|T2l has to be incorporated in the 
model as a constraint to select the physically relevant solution, as suggested in [DS96] and 
|LW07| . Nevertheless, the integrability of the fragmentation term is not a big deal since it 
can be obtained by imposing boundedness of a large enough moment of the initial data as it 
will be clear in the discussion below and as it appeared in [DS961 ILW07] . More interesting is 
how to interpret this in terms of boundary conditions; we shall discuss the point in Section 
|U (Note that in |LW07j the problem is completed with the boundary condition U(t, Xo) = 
while xo > 0, t(xq) > 0.) According to [DS96, LW07] we adopt the following definition. 

Definition 1 We say that the pair (U, V) is a "monomer preserving weak solution of the 
prion proliferation equations" with initial data (Uq,Vq) if it satisfies ([2]) and if for any 
<p G C%°((xq, 00)), we have 

poo 

U(t, x)(p(x) dx — / Uo(x)ip(x)dx 
Jo 

t roo rt /"OO 

/ n(x)U{s, x)<p(x) dx ds — / / f3(x)U(s, x)tp(x) dx ds 
Jo Jo Jo 

t roc />t /"oo ry 

+ 1 V(s) T(x)U(s,x)d x ip(x)dxds + 2 / f3(y)U(s,y) k(x, y)<p(x) dxdy ds, 

Jo Jo J x Q J x 

(13) 



and 



V{t)+ I xll(t,x)dx = V + xU (x)dx 



Xt — I jV(s) ds — / / xfi(x)U(s, x) dx ds 



(14) 



A break is necessary to discuss the functional framework to be used in Definition [T] We 
start with a set up of a few notation. We denote by M 1 (X) the set of bounded Radon 
measures on a borelian set X C R; M.\{X) stands for the positive cone in A4 1 (X). The 
space A'( 1 (X)_identifies as the dual of the space Cq(X) of continuous functions vanishing at 
infinity in A|^| endowed with the supremum norm, see [Mal82j . Given an interval / C R, 
we consider measure valued functions W : y G / 1— > W(y) G M. 1 (X). Denoting W(y,x) = 
W(y)(x), we say that W G C(I;M 1 (X) — weak — *), if, for any ip G Cq(X), the function 
U ^ f x ( f(x)W(y, x) dx is continuous on I. We are thus led to assume 

U G C([Q,T];M\ ([0,oo)) - weak-*), V G C([0,T]), 



1 4> £ Co(X) means hat is continuous and for any r) > 0, there exists a compact set K C X such that 
sup^^ |</>(a;)| < rj- We denote C C (X) the space of continuous functions with compact support in X. 
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with furthermore 

/>oc 

supp(J7(t, .)) C [xo, oo), / xU(t,x) dx < oo, 

which corresponds to the physical meaning of the unknowns. Hence, formula (fT3)l makes 
sense for continuous coefficients 

H, f3, t G CQx Q ,oo)). 
Concerning the fragmentation kernel, it suffices to suppose 

y h-> k(-,y) G C([x , oo); A^ + ([0, oo)) — weak — *). 

2 Main result 

2.1 Notations and rescaled equations 

We first rewrite system (p} in a dimensionless form, as done for instance in [CGPV02] (see 
also |LM07| ). We summarize here all the absolute constants that we will need in the sequel: 

• T characteristic time, 

• hi characteristic value for the concentration of polymers Ui, 

• V characteristic value for the concentration of monomers v, 

• T characteristic value for the polymerisation rate Tj, 

• B characteristic value for the fragmentation frequency Pi, 

• do characteristic value for the degradation frequency of polymers fii , 

• r characteristic value for the degradation frequency of monomers 7, 

• L characteristic value for the source term A, 

The dimensionless quantities are defined by 

- t v(tT) _ „ uAiT) - fa _ Ti 

t = j, , v{t) = — — , Ui{t) = - , Pi = —, n = — , 

Mi t ^ - 7 

tM = d- > A = Z' 7 = f 

We remind that kij is already dimensionless. The following dimensionless parameters appear 

LT 

a = — — , b = BT, c = TT, d = doT 

J (15) 
s=-, is = TTV. 

Omitting the overlines, the equation becomes 

— = aX — c^v — vsv + 2bs ^ ikijPjUj, 

j>n i<no (16) 

= -d^Ui - 6/3jUj - vv(T t Ui - n^Ui^x) + 2bY lj>i Pjh,jUj, for i > n . 
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The definition (fit)]) of the total mass in dimcnsionless form becomes 



v + s ^2 im = p. (17) 

The rationale motivating the scaling can be explained as follows. Let < e <C 1 be 
a parameter intended to tend to 0. We pass from the discrete model to the continuous 
model by associating to the u^s a stepwise constant function, constant on each interval 
(ei,e(i + 1)). Then sums over the index i will be interpreted as Ricmann sums which are 
expected to tend to integrals in the continuum limit while finite differences will give rise 
to derivatives. Having in mind the case of homogeneous division and polymerization rates 
/3(x) — x a , t(x) — x 6 , fi(x) = x m , which generalizes the constant-coefficient case proposed 
by [GPMW06J, and their discrete analogue /3j = i a , Ti = i e , /i, = i m , we shall assume that 
the rescaled coefficients f3i,/J.i,Ti fulfill ©. Therefore, we are led to set 

s = e 2 , 

so that (fT7|) becomes 

oo 

v + e ei Ui = p, (18) 

i=n 

to be compared to the definition of g in (|12|) . This scaling means that the typical concen- 
tration of any aggregate with size i > no is small compared to the monomers concentration, 
but the total mass of the aggregates is in the order of the mass of monomers. Next, we set 

a = l, b = e a , c = l, d = e m , v = e e ~ 1 . 

The rescaled equations read 

= A - >yv - e e+1 v Y, nui + 2e 2+a ^ ]T j i 

i> n o 0<no 

= -e m fj,iUi - e a (3 t Ui - e e ~ 1 v(T i u i - r,_iUj_i) + 2e a 2J Pjhjuj, for i > n . 

j>i 

(19) 

Eventually, the threshold value uq also depends on the scaling parameter and we assume 

limen (e) = a; > 0. (20) 



This choice is discussed in Section 15.31 

Equation (|19|) is completed by an initial data (u® ,£ ,v 0,£ ) verifying, for some constants 
Mo, p° , Mi +ty independent of e : 

oo 

v 0,e + £ 2 iu"-" = p° < +00, 

i=n (e) 

oo 

e Yl u ^' £ ^ A/ ° < +°°' (21) 



i=n (e) 

oo 



e 2+a Y i 1+<7 u° i ' £ < M 1+a < +oo, 1 + cr > max(l,a, 1 + m, 1 + 6). 

i=n (e) 

For any < T < oo, Theorem Q] guarantees the existence of a solution (uf,v e ) of (fT9|) . Let 
us set 

Xl(x) = X[ie,(i+l)e) (x) 
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with xa the indicator function of a set A. We introduce the piecewise constant function 

oo 

u E (t,x):= Y, «"(*)x?(aO 

i=n (e) 

On the same token, we associate the following functions to the coefficients 

oc oo , 

k £ (x,y) :=X)E-rX?(s)x?(tf), 
i=0 i=0 

oo oo oo 

//(*):= ]T £ m /iiXf(a:). TO ~ E ^ftxfW. ^(as) := £ e fl TiXf(z)- 

i=n (e) i=n (e) i=n (e) 

This choice is made so that for all y, k £ (-, y) is a probability measure on [0, j/]. 



2.2 Compactness assumptions on the coefficients 

For technical purposes we need further assumptions on the discrete coefficients. Let us 
collect them as follows: 

There exists K > such that 
f3i+i-/3i\ <Ki a ~ x 

n+i -n\ < Ki e ~\ 

For the fragmentation kernel we assume 



(22) 



where the exponents a, 9, m arc defined in i 
furthermore 

There exists K > such that for any i, j 

i-l p—1 i-1 p— 1 



< 



< K. 



(23) 



E E k r,3 + l E E fcr ^ 

p— r— p— r— 

These assumptions will be helpful for investigating the behavior of (p~9|) as e goes to 
since they provide compactness properties. We summarize these properties in the following 
lemmata. 



Lemma 1 Let (zt 



be a sequence of nonnegative real numbers verifying 



< zt < Ki K , 



\Zi + i - zA < Ki h 



for some K > and k > 0. For x > 0, we set z E (x) = £ K z iX[ei.e(i+i))( x ) ■ Then there 
exist a subsequence e n — > 0, and a continuous function z : x £ [0, oo) t— > z(x) such that z £n 
converges to z uniformly on [r,R] for any < r < R < oo. If k > 0, the convergence holds 
on [0, R] for any < R < oo and we have z(0) = 0. 

We shall apply this statement to the sequences f3 £ , fi £ ,T £ . A similar compactness property 
can be obtained for the fragmentation coefficients. 



Lemma 2 Let the coefficients ki t j satisfy Assumptions {Sp, (OJ) and Ii23\) . Then there ex- 
ist a subsequence (£n) JlgN o/nd k : y £ [0, oo) i— ► k(-,y) £ A4]j_([0, oo)) which belongs to 
C([0, oo); A4+([0, oo)) — weak — *) satisfying also f3|) and 0) (in their continuous version) 
and such that k £ " converges to k in the following sense: for every compactly supported 
smooth function if £ C%°([xq, oo)), denoting 



n (e n )e 



k £n (x,y)(p(x)dx, <f)(y) = k(x, y)<p(x) dx 



(24) 



we have (jf n — > <p uniformly locally in [xq,+oo). 

The detailed proofs of Lemma [T] and Lemma [2] are postponed to Appendix IA.2I 
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2.3 Main results 

Wc are now ready to state the main results of this article. 

Theorem 2 Assume (|5J) and (|22p . Suppose the fragmentation coefficient fulfill ([I])-© and 
(|23p . Then, there exist a subsequence, denoted (£n) ngN , continuous functions /j,,T,f3, and a 
nonnegative measure-valued function k verifying 0) and such that 

//" , t e " , f3 E " , k e " -> fj,, t, /3, k 

in the sense of Lemma[l] and Lemma\^ 

Let the initial data satisfy (|21[) . Then we can choose the subsequence (£ n ) nGN such that 
there exists (U, V) for which 

u E » ->> U, mC([0,T];7W 1 ([0,oo)) - weak-*)), 
ti E ™ — * V uniformly on [0,T]. 

We have xU(t,x) e ^([0,00)), the measure U(t,.) has its support included in [xo,+oo) 
for all time t>0, and (U, V) satisfies ([T3 ]l -([13 ]) . 

Theorem 3 The limit (U, V) exhibited in Theorem\^is a monomer preserving weak solution 
(i.e. satisfies also Equation in the following situations: 

i) Xq = and either 9 > (so that the limit r satisfies t(0) = 0), or the rates Ti = r are 
constant. 

ii) Xq > and the discrete fragmentation coefficients fulfill the following strengthened 
assumption: for any i,j we have 



i'j + l ki> j 



< j, h,j < j. (25) 



3 Moment estimates 

We start by establishing a priori estimates uniformly with respect to e. These estimates 
will induce compactness properties on the sequence of solutions. As described in [Lau02j 
for general coagulation fragmentation models, the model has the property of propagating 
moments. 

Lemma 3 Let the assumptions of Theorem^ be fulfilled. Then for any T > 0. there exists 
a constant C < oo which only depends on M , M 1+cr , K and T , such that for any e > 0: 

POO 

sup / (1 + x + x 1+a )u E (t,x)dx < C, 0<v E (t)<C 
te[o,r] Jo 

Proof. For r > 0, we denote 

oo 

M E (t)=eJ2(ie) r «!(*)• 

i=na 

As in [CGPV02] . we can notice that 

(^yu 6 (t,x)dx<M E (t)< ( x r u £ (t,x)dx. 
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Therefore, we only need to control Mq(£), M-f(t) and M[ +a (t). We notice the obvious but 
useful inequality, for < r < 1 + a. 

iie) r < l + (ie) 1+a , 

and therefore, 

|M r e | < \M$\ + \Mf +<7 \. 

In the sequel, we use alternatively two equivalent discrete weak formulations of Equation 
(TT9"]) in the spirit of |LW07| . We multiply the second equation of P^|l by <pi and summing 
over i, we first obtain 

, oo oo oo 

uftpi = -e m £ WiVi ~ £ " E PWiVi 



dt 

i=n i=n 
oo 



-e 9 - 1 v 6 ( r *< ~ n-iUi-ityi + 2e ° E fi^Pjki 

i— no i— no j>i 

oo oo oo 

= -s m £ \nu\ipi - e a £ Piulipi + e 6 - 1 £ nv%(<p i+1 - <p 

i—rtQ i—uq i—riQ 

oo 

+2e a ^X '''■>'''•■•",■ 

i=n j>i 

Using the properties of kij, we rewrite the fragmentation terms as follows 



(26) 



j=n + l i=l ] 

oo j—1 oo no — 1 

2 X E ik '-> j J U J 2 J 1 2 E E ***,J& u ?y +0n o < o <Pno, 

j=n a + l i=n j=n +l i=l 

i-i 



2 X] ^E^Aj'^ = 2 E ^2 ik idPj u j-f- 

i=n j>i j=n Q + li=n 

By using JZ]), we obtain 

oo oo oo no — 1 

2 y. v»S - £ = - 2 E E ifej&^y 



i=n j>i i=n j="o i=l 

oo 2-1 



2 E E*^ 

j=n +l i=n 



i 3 



Replacing in the weak formulation we get 



■jt E u ^ = ~ e ™ E M^f^ + e<? 1?;e E TiU ^ ~ ^) 

i=no i—UQ i—riQ 

oo j—1 / \ oo n - 1 

j=n + l i=ri ^ ' j=n a i=l 

(27) 

This last formulation makes the estimates straightforward (the computations are formal but 
can be understood as uniform bounds on solutions of truncated systems and therefore on 
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any admissible solution). Taking <pi = is, we obtain the first moment, that is, the previously 
seen mass balance: 

— ( v E + e 2 iu i) = "7"" ~ £2+m Y ^< + A < A. (28) 



i— uq ' %—riQ 



Therefore, we get (uf and v e are nonnegative) 

< v £ (t) + M[(t) <p° + XT for < t < T < 00 

and 

/ e 2+m Mi ^( S) x ) d s < p° + XT for < t < T < 00. 
Jo 

i=n 

To obtain an estimate on the Oth order moment, we take ifi = e. The term with r, vanishes. 
Considering only the nonnegative part of the derivative, we derive from (|27|) 

-M £ (t) < 2e 1+a Yl Y 

j=n Q +l i=n 
00 

< 2e 1+a ^ 2KM a(t)- 

j=na + l 

To give the bound on the (l+er)th moment, we choose ipi = e(ei) 1+a in the weak formulation. 
Thanks to the mean value inequality, we have 

((e(* + 1)) 1+CT - (ei) 1+a ) < (1 + a)e{e(i + 1)) CT < (1 + a)2 (r e(ei) (r , 

therefore (f27|) yields 

j 00 00 

-Mf +CT (0 + e 1+m Y ^(£*) 1+ X < + ^ £ e^ufe^)", 

i=no i—riQ 

< K{p° + \T){l + a)2°Ml +a {t). 
Since < 6» < 1, denoting C = max (A' (p° + AT)(1 + o-)2 CT , 2A), it leads to 

_d 

df 



Af e (t) + Mf, Ji) < C M«(t) + M$ +S (t) <2C M £ {t) + Ml +a {t) , 



and we conclude by the Gronwall lemma. It ends the proof of Lemma [3] □ 

Hereafter, we denote by C a constant depending only on T, Mo, p°, Mi +cr , A and A such 
that 

M^,v s ,M?,M[ +a <C. 

Lemma 4 Under the assumptions of Lemma the sequence of monomers concentration 
( v£ )e>0 is equicontinuous on [0,T]. 

Proof. We use the estimates of Lemma [3] to evaluate the derivative of v £ . We recall the 
equation satisfied by v £ 

^ = A - iv £ + e 1+ V £ nu! + 2e 2 +« £ Y jfc^f, 

!>n j<no 
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which implies 



dv £ 



dt 



< A + 7 C + KC 2 + 2sn Q {e) KM e a . 



Since the sequence (M^) £>0 is uniformly bounded with respect to e by Lemma [3] (recall 
that a < 1 + a), the sequence (i> £ ) £>Q satisfies a uniform Lipschitz criterion on [0, T]. This 
concludes the proof of Lemma 01 □ 

Proof of Theorem [2j By the Arzela-Ascoli theorem and Lemma [4] there exists a 
function V G C([0, T]) and a subsequence that we still denote v £ such that 

v £ (t) — >V(t) in C([0,T]). 

In the same way, the moment estimates of Lemma [3] give uniform boundedness for (1 + x + 
.t 1+ct )w £ in A^dfLoo)). Pick a function ip G C~([0,oo)). We define 



so 



that V u?93f = / 
Jo 



(i+l)e 

9?(x) da;, 

it e (i, x)i^(x) da;, and also for y G [j£, (j + l)e[ 



fc 6 (x,y)tp(x)dx — I k e {x,j£)(p{x)dx = 2_jkij — 



Thanks to the moment estimates of Lemma [3j and using (|26[) , we have 



(1 

~dl 



u £ (t, x)tp(x) dx 



<c(|M|oo + 11^1100) 



and 



u £ (t, x)ip(x) dx 



<C\\<p\\ 



for some constant C depending only on K, Mq, Mi +a , A, T. Therefore, for any function 
ip G C£°([0, 00)), the integral J u £ (-, x)p{x) dx is equibounded and cquicontinuous. Using a 
density argument, we can extend this property to tp G Co([0, 00)), the space of continuous 
functions on [0, 00) that tend to at infinity. This means that (/ °° u £ {., x)<p(x) dx) £ belongs 
to a compact set of C(0, T). As in |CGPV02j . by using the separability of Co ([0, 00)) and the 
Cantor diagonal process, we can extract a subsequence u £n and U G C([0, T\;M. ([0, 00)) — 
weak — *) , such that the following convergence 



u En (t, x)tp{x) dx 



U (t, x)p(x) dx, 



as e n — > 0, holds uniformly on [0, T], for any 93 G Co([0, 00)). As u e (i, x) — for a; < eno(e), 
we check that U(t, .) has its support in [xq, oof. It remains to prove that ({/, V) satisfies (|13|) 

JH. 

Let <p be a smooth function supported in [5, M] with a:o < 5 < M < +00, choosing 
eno(e) + 2e < <5 (what is possible due to (HO)))- By using Lemma [1] and Lemma we check 
that, for a suitable subsequence, one has 



fi 6 " (x)u £n (t, x)cp(x) dx 
/3 £r> (x)u £ ™ (t, x)<p(x) dx 
r £n (x)u £n (t, x)ip(x) dx 



fj,(x)U(t, x)(p(x) dx, 
f3(x)U (t, x)ip{x) dx, 
t{x)U (t, x)p(x) dx, 



(29) 
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uniformly on [0, T\. Equation (|26|) can be recast in the following integral form 

J />oc />oo />OC 

— / u £ (t, x)ip(x) dx = — / [i £ u £ (t, x)tp(x) dx — / T £ u £ A £ ip(x) dx 

Jo Jxo JO 

/>oo poo />oo 

- (3 £ u £ {t,x) V {x)dx + 2 / <p(x)p E (y)u e (t,y)k e (x,y)dxdy 

JO Jo Jx 

where we have defined 

A>(x) = j (t+l)£ ^ S + £) '^ (S) da, for x e [fe, (i + 1)4 



(30) 



The first and third terms are treated in (f2U)) . Using (|2T))) again and remarking that 
|A%t) - ip'(x)\ < ellv/'Hoo, we have 

r e "(x)u 6 "(t,x)A e "< y 5(a;)da; ► / T(x)U(t,x)(p'(x) dx, (31) 

uniformly on [0, T]. Let us now study the convergence of the last term in (|30"|) . To this end, 
we use the notation and <p £ as defined in |24|) of Lemma [2] and we rewrite 

oo />y />oo 

/ ^(x)k £ {x,y)u £ {t 1 y)f3 £ {y)dxdy = 2 u e {t,y)l3 e {y)(jf{y)dy. 

X J Xq Jx 

Owing to (f2"3")) we use Lemma [5] which leads to 

<jf n ^> <j) uniformly on [xo,-R] for any R > 0, 

and thus also 

/3 E "f" > @<f) uniformly on [xo,-R] for any R > 0, 

for a suitable subsequence. Finally, we observe that e " and therefore 4> are bounded by 
H^Hoo. Thus, by using the boundedncss of the higher order moments of u £ in Lcmma[3]with 
1 + a > a, we show that the fragmentation term passes to the limit (see Lemma [5] in the 
Appendix). We finally arrive at 

/•OO 

U(t, x)<p(x) dx — / f/(0, x)ip(x) dx 

J x 

t />oo ft />oo 

/ fiU(t, x)ip(x) dx — / V(s) / t(x)U(s, x)ip' (x) dx 

Jx JO Jx 

t />oo r-t />0O f-y 

\ 0(x)U(s,x)(t,x)<p(x)dx + 2 / p{y)U(s,y) (p(x)k(x,y) dxdy, 

'0 Jx Jo Jx Jo 

(32) 

which is the weak formulation (fT3l). Moreover, (l28l) recasts as 



v £ (t)+ / e £ (x)u £ (t,x)dx = v°' e (t)+ / e £ (x)u £ (0,x)dx 



+\t — 7 / v £ (s)ds- / / e e (x) / u e (x)M £ (s,x) dxds 
Jo Jo 
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where 

oo 

e £ (x) = Y ei X[ie,(i+1) 



i=0 

Clearly e £ {x) converges to x uniformly. Using the moment estimate in Lemma[31 with a > 0, 
we obtain 

poo poo 

v £n (t) + e £n (x)u £ "(t,x)dx >V(t)+ xll(t,x)dx 

Jo Jo 

uniformly on [0, T] as well as 

t pOO r-t />0O 

e £n (x)/j, £rl (x)u £ " (s,x) dxds > / / xfi(x)U(s, x) dxds. 

Jo Jo 

(We refer again to Lemma or a slight adaptation of it.) As e n — > we are thus led to 
©• □ 

Proof of Theorem [3j We rewrite the rescaled ODE as 

dv e f°° r°° r° 

— = \- lv e - T £ (x)u £ (t,x)dx + 2 p £ (y)u £ {t,y) e e (x)k £ {x,y)dx, 

at Jn a e Jn e JO 

Depending on the value of xq, we have to care about the last term (xq > 0) or the next two 
last term (xq = 0). As already remarked in the proof of Lemma 01 in case where xq = 0, 
the fragmentation term can be dominated by 



2e2+Q 3kj,iPiU e < 2eno{e) KM £ a . 



i>no 3< n 

Hence this contribution vanishes as e goes to when lim e ^o £no(s) = xq = 0. Nevertheless 
for case i) we still have to justify that J °° T £ (x)u £ (t, x) dx passes to the limit. We get 

/>oo 

T £n u En (t,x)dx > rU(t,x)dx, in C([0,T}) (33) 

"o(Eb)e„ En^O J Xo 

by using the strengthened assumption < 9 < 1 in (f2"Tj) . Indeed it implies that t £ (x) 
converges uniformly to t(x) on any compact set [0, R] while these functions do not grow 
faster than x at infinity. We can thus use Lemma [5] to conclude. 

In the situation ii). another difficulty comes from the fragmentation term since we have 
to prove that 

n (e n )e n /■oo rx 

e £ "(x)k £ »(x,y)p £ »(y)u £ "(t,y)dxdy 2 / / xk(x, y)j3(y)U(t, y) dx dy. 

n (E„)e„ JO e ™ _> " Jx a JO 

The stronger compactness assumptions (|2"5)l are basically Ascoli-type assumptions on the 
repartition function associated to the kernels k £ . Denoting, in a similar manner to Appendix 
1X21 

pX pX 

F £ (x,y)= k £ (z,y)dz, G £ (x,y) = F £ (z,y)dz, 



Lemma [TU] (see Appendix IA.2[) ensures that F £ — > F uniformly on compact sets of R+ x 
[xo, +oo). We also make the remark that 

rn e rnoe 

/ e £ (x)k £ (x, y) dx — / xk £ (x,y)dx < e, 
Jo Jo 
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r n e 



10 



x—noe 



x=0 



"of 



xk £ {x,y)dx = xF £ (x,y) F £ (x , y) = (n e) F £ (n e , y) - G £ (n e , y) . 



o 



Thanks to Lemma 1101 we know that the concerned quantities are uniformly bounded and 
converge uniformly on compact sets, so that 

/ e £ " (x)k £ " (x,y) dx ► / xk(x,y)dx uniformly on compact sets. 

Jo e™^0 J Q 

And as before this is sufficient to prove that 

poo pn e n roc rx 

2 l (3 e "(y)u £n (t,y) e £n {x)k £n {x 1 y)dxdy -» 2 / 0(y)u(t,y) xk(x,y)dxdy. 

Jn e n JO £n^0 J XQ J 

□ 

4 Boundary Condition for the Continuous System 

The discrete system (TTJ) only needs an initial condition prescribing the u^s and v at time 
t = to be well-posed, as Theorem [T] states. It is different for the continuous system ©(O: 
a boundary condition at x = xq is needed when t(xq) > (in which case the characteristics 
associated to the "velocity" r are "incoming"). Even when t(xq) = 0, difficulties might 
arise when x i— > t(x) is not regular enough to define the associated characteristics. However, 
according to the analysis of [DS96, LW07], we have seen that the notion of "monomer 
preserving solution" appears naturally, inserting (|12[) as a constraint. It leads to the question 
of deciding how this condition is related to J2]) and (|3|) and to determine the corresponding 
boundary condition to be used at x = xq. 

Let (U, V) be a "monomer preserving" solution. In this section we do not care about reg- 
ularity requirement, and we perform several manipulations on the solution (that is assuming 
all the necessary integrability conditions). We suppose that the kernel k splits into a Dirac 
mass at x = xq and a measure which is diffuse at xq: 

k(x,y) = l(x,y) + 6(x = x+)ip + {y) + 5(x = x^)tp~(y), 

where for any y > 0, f^°_ n l{x, y) dx — > as r\ goes to 0. We have to distinguish between x^ 
and Xq since their biological and mathematical interpretation is different: the Dirac mass at 
Xq means that polymers of size xq are formed, whereas the Dirac mass at Xq is interpreted 
as breakages of polymers of size Xq going back to the monomers compartment V. As shown 
below, the mathematical treatment of each is different. The time derivative of (p~4|) leads to 

d dV /"OO Q l>OC 

— g=— — h / x— U(t, x) dx = — / xfi(x)U(t, x) dx + A — 
dt dt J XQ dt J XQ 

In the left hand side, we can compute the derivative of the moment of U by using ©. We 
get 

^ /"OO POO f-OQ 

— I xU(t, x)dx = — I x(3U dx — xpJJ dx 

d Jx ° Jx ° r oo a Jxo ,00 

-V x—(rU)dx + 2 xl l(x,y)p(y)U(t,y)dydx. 

J Xq ®% J Xq Jx 

In this equation, since ([3]) is only written for x > xq, only the diffuse part of the kernel k 
appears. Integrating by parts, the convection term yields 

/ x— — (tU) dx = — xoT(xo)U(t, xq) — / rlldx. 

Jxn V X Jxn 



pOG /"GO 
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Now wc use (|2J), which writes 

JV pOO pOO pXg 

— =X-jV-V r{x)U(t,x)dx +2 / y l(y,x)P{x)U(t,x)dydx 

dt J Xq Jx=x a J y=0 

/>oo 

+2x Q / i/)~(x)/3(x)U(t,x)dx 

J X — Xq 

then we obtain 

x V(t)T(x )U(t,x )- x(3(x)U(t,x) dx+2 f3(x)U(t,x)( yl(y, x) dy + x ip~{x) ) da; = 0. 
However, ([7]) is interpreted as 

2 / vKy,x)dy J r2x G X[x ,oo){x)'tp~{x) + 2x(iX{x a ,oo){x)%l) + (x) = x. 
Jo 

We are thus led to the relation: 

.t (v(t)T(x )U(t, x ) - 2 J i)+{x)P(x)U(t, x) dx^j = 
which suggests the boundary condition 

/>oo 

x VT(x )U(t,x ) =2x / V> + (^)/?0E)lT(M)d2;. (34) 



(Note that written in this way it makes sense also when xq = 0.) 

When xq > 0, the above calculation gives solid intuitive ground to choose Equation (|34|) 
as a boundary condition, defining the incoming flux by means of a weighted average of the 
solution over the size variable. In particular if the Dirac part vanishes we obtain 

VT(x o )U(t,x o )=0, 

the boundary condition proposed in |GPMW06] . for constant coefficient r. It is also the 
boundary condition used in |LW07j . 

If xo = 0, the problem is still harder, since Equation (|34j) is empty. Dividing it by xq > 
and passing formally to the limit would however give: 

poo 

VY(0)t/(0) = 2 / iP+(x)P(x)U{t,x)dx. (35) 



Here again, it generalizes what has been proposed in GPMW06 for t constant and k(x,y) = 
yXx<y, though without any rigorous justification, and if ip + = it imposes a vanishing 
incoming flux. 

5 Discussion on the parameters and choice for e 
5.1 Orders of magnitude 

A biological discussion upon the parameter values can be found in [Lcn09j and is based on 
[MJN991 IMGA05| and references therein. 

To carry out the previous scaling limit theorem, we made the following assumptions: 

s = — = e 2 , v = -, lim £ _ > o£"o(e) = %0i rj = a = c = d=\. 
V e 
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Let us denote io the average size of polymers. Even if there still exists much uncertainty 
upon its value, we can estimate that the typical size of polymers ranges between 15 and 
1000, so we can write 

1 

£l = — < 1. 

to 

It is also known that the "conversion rate" of PrPc is around 5 to 10% at most (depending on 
the experiment, on the stage of the disease, etc) ; it means that the mass of proteins present 
in the monomeric form is much larger than the mass of proteins involved in polymers. In 
terms of characteristic values, it writes 



e 2 = — « 1. 




Finally, we get: 



Hence, it legitimates the assumption on the parameters s and e. Concerning the parameter 
a, we have a = ^ ~ , which is in the order of 1. We have only do < 5.10 -2 : this should 
lead to neglect the degradation rate of polymers and simplify the equation. 

For the fragmentation frequency, it is in the order of the exponential growth rate, found 
experimentally to be in the order of 0.1; in the case of Masel's articles jMJ N99. MGA05], it is 
supposed that a = 1, so it seems relevant (it leads to a fragmentation frequency in the order 
of e). However, it has to be precisely compared to the other small parameters which are 
given by the typical size i$ and the conversion rate to justify the approximation. Moreover, 
the assumption of a linear fragmentation kernel (3 has to be confrounted to experiments. 

Concerning the aggregation rate T, and its related parameter v = tV, as shown in [Len09j , 
in most cases we have ^ in the range of [0.01, 0.1], so it seems justified to suppose it small ; 
what has to be explored is its link with the other previously seen small parameters. 

To conclude (or open the debate), it seems that each specific experiment, like PMCA 
protocole, in vitro or in vivo measures, or yet for the case of recombinant PrP (see |Rez08] ) . 
the orders of magnitude of each parameter should be carefully estimated, in order to adapt 
the previous model and stick to the biological reality - which proves to be very different in 
in vivo, ex vivo or in vitro situations, or yet at the beginning (when there are still very few 
polymers) and at the end of experiences. The following discussion illustrates this idea, and 
gives some possible extensions to the previously seen models. 



5.2 Discussion on the fragmentation rates kij 

To illustrate the central importance of a good estimate of the orders of magnitude, we 
exhibit here a case where the limit is not the continuous System ©CI]), but another one. 
Our calculation is formal, but a complete proof can be deduced from what preceeds and 
from |CGPV02j . 

Let us take, instead of b = e a : 

b = e a -\ 

and suppose also that the fragmentation kernel verifies: 

h,i = h-i,i = ^(l - en), hj = . < i < j - 2. 

It means that the polymers are much more likely to break at their ends than in the middle of 
their chain. In that case, under Assumption ([23)) on k®j and (|22|) on tj and j3j, if a— 1 < 1+a, 



1G 



the limit system writes: 



dv 

= X — "/v + v I T(x)U(t,x)dx 



/3(x)U(t, x) dx + 2 / / yk(y,x)r(x)U(t,x)dydx, 



dt 

du d d f°° 

— = -tx(x)U(t,x)-r(x)U(t,x)-v— {tU) + —((3u) + 2 I r(y)k(x,y)U(t,y) dy. 

(36) 

Notice also that System ([2]) ([3]) includes the case of "renewal" type equations (refer to |Per07j 
for instance), meaning that the ends of the polymers are more likely to break. For instance, 
if we have, in the above setting: 



2 J 

then Equation ([36]) remains valid, but we have to write the boundary condition as: 

t(x = 0)U(x = 0) = J m(y)U(t, y) dy. (37) 

This indeed can also be written as the measure of by d[i y = k(x, y) dx : m(y) = /ij,({0}). 

Both of these cases mean that the ends of polymers are more likely to break. What 
changes is the order of magnitude of what we mean by "more likely to break": is it in 
the order of i, in which case System ©O is valid but with a boundary condition of type 
(f37|) ? Or is the difference of the order of i, in which case Equation (f36|) is more likely ? 
Refer to |Len09j for a more complete investigation of what model should be used in what 
experimental context. 

5.3 Discussion on the minimal size uq 

We have seen above that to have xo = 0, it suffices to make Assumption (j2T))) . Having also 
seen that the typical size io is large, and that 



iQm\V m\V 
it is in any case valid to suppose that 

— =e c , 0<c<2. 

Hence, Assumption (f2T)| can be reformulated as: 

n <4. (38) 

For c = 1, it means no <C «0j which is true. On the contrary, if we suppose that xq > 0, 

1 

it means that n « ig : in most cases, where for instance io = 100 or i a = 1000, it seems 
irrelevant. 
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A Appendix 

A.l Compactness of the coefficients 

Proof of Lemma [TJ We refer to [CGPV02j for the case k = 0. We prove here the case 
k > 0. First, we show that z e is close to a subsequence satisfying the requirements of the 
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Arzcla-Ascoli theorem on [r, R] . We define z e by 



We have 



z £ (x) = e K Zi + — (x - ie) for is < x < (i + l)e. 



\~z s (x)-z%x)\ = \s K Zi+1 Zi {x-ie)\, 



< £ K \z i+1 -Zi\, 

< eKiei)*- 1 < 2£(Kr K ~ 1 

Furthermore z e has a bounded derivative since 

dz E 



KR K 



dx 



Z l+ l - Zi 

= £ , 

£ 

< K{ei) K -\ 

< Kr*' 1 + KR K ~ 1 . 

Therefore, the family z e satisfies the requirements of Arzela Ascoli theorem for any interval 
[r, R] with < r < R < +00. We can extract a subsequence converging uniformly to z. The 
limit is continuous and satisfies z(x) < Kx K . When n > the convergence extends on [0, R] 
owing to the remark 

sup \{z e - z){x)\ < 2Kr. 

x£[0,r] 

This concludes the proof. □ 

During the proof of Theorem [2] and Theorem [3] we made repeated use of the following 
claim. 

Lemma 5 Let z n converge to a continuous function z uniformly on [0, M] for any < 
M < 00, with \z„(x)\ < K(l + x K ). Let (w n ) ngN be a sequence of integrable functions which 
converges to U weakly-* in Al 1 ([0,oo)). We suppose furthermore that 

p oc 

sup / (1 + x e )\u n (x)\ dx < C < 00. 



Assuming < k < I, we have 

z n (x)u n (x) dx 



z{x)U{x) dx. 



Proof. Let C G C c °°([0,oo)) such that < £(x) < 1, £(x) = 1 on [0, R], < R < 00 and 
supp(C) C [0,2R]. We split 



z n (x)u n (x) dx — / z(x)U{x)dx 

(z n {x)u n {x) - z(x)U(x)) (C(x) + 1 - C(x)) dx 

/•OO 

z„(x) — z(x)\ \u„(x)\ C(x)dx+ / z(x)C(x) (u n (x) - U(x)) dx 

Jo 

\z n (x)u n (x) - z(a:)i7(a;)| (l - ((x)) dx. 



< 



The last integral can be dominated by 
'l + y K ' 



K sup 



v >r V 1 + y 



sup/ (1 + a:')(|u„(a:)| + |I7(a:)|) dx 
n Jo 
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Since < k < £, this contribution can be made arbitrarily small by choosing R large enough, 
uniformly with respect to n. Moreover, we clearly have 

/>oo />oo 

0< / \z n (x) - z(x)\ \u n (x)\ ((x)dx < sup \z n (x) - z(x)\ sup/ \u m (x)\dx >0 

JO 0<x<2R m Jo jwoo 

and of course 

/■oo 

/ z(x)C(x) (u n {x) — u(x)) dx > 0. 

Jo ™^°° 

Combining all together these informations ends the proof. □ 
A. 2 Compactness of the fragmentation kernel 

We look for conditions on the coefficients guaranteeing some compactness of k £ . We use a 
few classical results of convergence of probability measures (see |Bil99j for instance). Let us 
introduce a few notations. Given a probability-measure-valued function y £ R *—>■ k(.,y) € 
we denote F(.,y) its repartition function: F(x, y) = J_ k(s, y) ds and G(x, y) the 
function F(z,y)dz. We shall deduce the compactness of k e from the compactness of 
the associated G £ . To this end, we need several elementary statements. 

Lemma 6 Let {P™, n € N} be a family of probability measures on R ; having their support 
included in some interval [a,b]. We denote F n the repartition function of P n , and G n the 
functions defined by J_ F n (s) ds. The following assertions are equivalent: 

1. P n P weakly (i.e., V / G C b (M), P n f Pf) 

2. F n (x) — > F(x) for all x at which F is continuous 

3. G n — > G uniformly locally 

Lemma 7 (Conditions for F) Let F be a nondecreasing function on R. There exists a 
unique probability measure P on R, such that F(x) = P(] — oo, a;]), iff 

• F is rightcontinuous everywhere, 

• lim^-oo F(x) = 0, lim^^+oo F(x) = 1. 

Furthermore P has its support included in [a, b] iff F = on] — oo, a[ and -F(6) = 1. 

Lemma 8 (Conditions for G) LetG be a convex function onR. There exists a probability 
measure P on R, having its support included in [a, b], such that G{x) — J_ cx> F(s) ds, where 
F(x)=P(]-o ,x}), iff 

• G is increasing, 

• forx> b, G(x) = G(b) + x, 

• G = on ] — oo, a]. 

Corollary 1 Let (G n ) a sequence satisfying the assumptions of lemma [3 Suppose 
G n — ► G uniformly locally on R, f/ien G afeo satisfy these assumptions and we have P n — > P 
weakly. 

Proof. We define the function F as F(x) = lim^^o+ G ( 2 + a )~ G '( a:: ) ; ^ j s then easy to check 
that F satisfies assumptions of lemma [3 and G(x) = F(s) ds. □ 

Proof of Lemma [2j We prove the following result, which contains Lemma O 
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Lemma 9 Suppose that the discrete coefficients satisfy 123\) . Then there exist a subsequence 
e n and k € C([0, oo), A^+([0, oo)) — weak — *) such that 

• k satisfies (6|/, {2P (and therefore ffi)), 

• for every y > 0, fc £ ™(., y) — ► k(.,y) in law, 

• for every ip € C£°([0, oo)), <j) E " — > <f> uniformly on [0, R] for any < R < oo. 

For any y > 0, fc e (a;, y) da; defines a probability measure on [0, oo), supported in [0, y]. We 
set F e (x,y) = f*k E (z,y)dz and G 6 (x,y) = f*F e (z,y)dz. Let tp € C^°(R* + ). We start by 
rewriting, owing to integration by parts, 

rv ry 
4> e {y) = <p(y) - / F 6n (x,y)<fi'(x)dx = (p(y) - G E (y,y)(p'(y) + / G E (x, y)<p"(x) dx, 
Jo Jo 

where we used the fact that F e (y,y) = f Q k e (z,y)dz = 1. The proof is based on the 
following argument: G £ is close to a G e which satisfies the assumptions of the Arzela-Ascoli 
theorem. Given x,y > and e > 0, i,j denote the integers satisfying x £ + l)e[, 

2/ G + an d a short computation leads to 

F e (x,je) = S it j + — — ki tj , 

G e (x, je) = e^5 P)j + (s - + ^S itj + ~ £ ' 

p=0 



where 



We define 



and we have 



Observe that 



Si t j — 22 k 



i-l 



r=0 



~k s (x,y) = (j + 1)£ V k s (xje) + V -^k^x, (j + l)e) 

e e 



&{x,y) = (J + 1)£ V G*{x,je) + V -^G*{x, (j + l)e). 

e e 



\G s (x,y)-G%x,y)\ = IL^f^, (j + l)e) - G s (x, je)\ 

£ 

< e^(Sp,j + i - S pd ) + (x- ie)(Si, j+1 - S itj ) + |(Si,j+i - S^) 

Due to (O, we have < fcjj < 1 and thus \hj + i — kij\ < 1. Similarly < 5jj < 1 and 
|Sij+i — Si.j \ < 1. Hence, since (|2"5j) can also be written 

i-l 
p=0 

it allows us to obtain 

\G s (x, y) - G £ (x, y)\ < e{K + 1 + 1/2 + 1/2) = e(K + 2). 
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We also deduce that 



\d v G e {x,y)\ 



while 

Moreover, we have 



\G s {x,je)-G e {x,{j + l)e)\ 



\d x &(x,y)\<l. 



\&(x,y)\<2e(t + 2) 



<K + 2 



which is bounded uniformly with respect to e and < x,y < R < oo. As a consequence of 
the Arzela-Ascoli theorem we deduce that, for a subsequence, G £n converges uniformly to a 
continuous function G(x,y) on [0,R] x [0,R] for any < R < oo. It follows that 

e »(j/) xj)(y) = (p(y)-G(y,y)ip'(y)+ [ G(x,y)ip"(x)dx 



uniformly on [0,i?]. We conclude by applying Lemma [8] to the function x G(x,y), with 
y > fixed. □ 



Lemma 10 Suppose that the discrete coefficients satisfy 1125)) . Then F e and G s are uni- 
formly bounded and converge (up to a subsequence) uniformly on compact sets. 



Proof. Assumption (f2l>)) rewrites 



< 



K 



K 



J 3 
so, with the same notation for F e as for G 6 and fc e , we get 

y - j'e 



\F £ (x,y)-F s (x, ] 



F e (x, (j + 1 e) - F E (x, je) < — < — e, 



together (the computations are very similar to those above) with 



and 



\d x F"\ < 



1 



K K 2K 



< — < 



< 



je uqe xq 



IK 

d v F^\ < -|F £ (x, (j + l)e) - F^x, je)\ < — 



e y 

which leads to Ascoli assumptions and therefore the suitable compactness. □ 

With such assumptions, we can take into account any k of the form k(x, y)dx = ^ko(x/y)dx, 
including Dirac mass. If we consider such a distribution on [0, 1] (taken symmetric), then 
tie ki^ 

i - 1 



we can define fcjj as 



ki j — k() 



i-i'i-i 



with Sj the Kronccker symbol. With these notations, we have for p > j — 2, 



( 



i=0 



0, - 



and Sj—ij = Sjj = 1, which leads to 



p p 
j ' j - 1 
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3— 1 

3 

as < i < j, we have for any p < i, 



p-1 p 



j - 1 i 

the intervals [j, [ and [^y^, j5t [ are disjoint. This leads to 

- s •»( u [?■ AD + W u{f})+ M u { A» « 

which gives the criterion (|23|) . The limit is then obviously given by k(x, y)dx = ^ko(x/y)dx. 
A. 3 Discrete system 

We discuss here briefly the existence theorem for the discrete system. It is mainly an 
adaptation of theorem 5.1 in |Lau02j . We define the truncated system. Let N > no, 
consider the system 

, JV-1 N 

— = A - JV - V ^ T i U * + 2 E E ik hjPj U 3> 
i=n j=n i<n 

du N 

= ^^no Un o ~ Pn u no ~ VT n + 2 ^ ] f3jk n(): jUj , for , 

= — mm - PiUi ~ v{nui - n-iUi-i) + 2 E Pjk hj Uj, for n a < i < N, 

j=i+l 

— 7— = -fJ-NUN - PNUN + UTJV-IUJV-I- 
di 

(39) 

Existence, uniqueness and nonnegativity are immediate, we have immediately the weak 
formulation 

. / N x N JV-1 

d 



(If 



«(*)^ + Ui< ^* ) = ^ ~ - v e + ^ e TiUi((p i+ i - if, - -0) 

i— no 2— no 

JV 3-1 , \ 

+2 £ e^^-(t-t) 



j=n +l i=n 
N n -l 

+ 2 E E iKiPjUity-j)- 

j=n i=l 

(40) 

Let us denote C/^ the infinite sequence of functions defined by = u" if no < i < 
jV, t//^ = otherwise. The weak formulation gives moment estimates (and the moment 
estimates done in section 4 can then be thought as uniform bounds on truncated systems). 
This model has the property of propagating moments. 

With this type of initial condition, the proof of existence is based on the Ascoli theorem 
for the continuous functions . Thanks all the moments controlled on the initial data and 
the nice property of propagation of moments, we have bounds on the derivative of v N , 
and therefore, we can extract convergent subsequence. The limit satisfies the equation in an 
integral form (see [BC90] for a definition). For proving uniqueness, the procedure exposed 
in |BC90[|Lau02| applies, with a small modification due to death rates (the condition on the 
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moment of order 1 + m for the initial data insures the convergence of ^2 ifiiUi). 

Acknowledgments. We thank Pavel Dubovski for kind comments and suggestions about 
this work, and Natacha Lenuzza for fruitful discussions and help about orders of magnitude 
of the parameters. Part of this work has been supported by the ANR grant TOPPAZ. 



References 

[BC90] J. M. Ball and J. Carr. The discrete coagulation-fragmentation equations: exis- 
tence, uniqueness, and density conservation. J. Statist. Phys., 61(l-2):203-234, 
1990. 

[BCP86] J. M. Ball, J. Carr, and O. Penrose. The Becker-Doring cluster equations: 
basic properties and asymptotic behaviour of solutions. Comm. Math. Phys., 
104(4):657-692, 1986. 

[Bil99] P. Billingsley. Convergence of probability measures. Wiley Series in Probability 
and Statistics: Probability and Statistics. John Wiley & Sons Inc., New York, 
second edition, 1999. 

[CGPV02] J.-F. Collet, T. Goudon, F. Poupaud, and A. Vasseur. The Becker-Doring 
system and its Lifshitz-Slyozov limit. SIAM J. on Appl. Math., 62(5):1488- 
1500, 2002. 

[CLD+08] V Calvez, N Lenuzza, M Doumic, J-P Deslys, F Mouthon, and B Perthame. 

Prion dynamic with size dependency - strain phenomena. J. of Biol. Dyn., in 
press, 2008. 

[CLO+09] V. Calvez, V. Lenuzza, D. Oclz, J.-P. Deslys, P. Laurent, F. Mouthon, and 
B. Perthame. Size distribution dependence of prion aggregates infectivity. Math. 
Biosci, 1:88-99, 2009. 

[DC95] F.P. Da Costa. Existence and uniqueness of density conserving solutions to the 
coagulation-fragmentation equations with strong fragmentation. J. of Math. 
Anal, and Appl, 192:892-914, 1995. 

[DS96] P. Dubovskii and I. Stewart. Existence, uniqueness and mass-conservation for 
the coagulation- fragmentation equations. Math. Methods in the Appl. Sci., 
19:571-591, 1996. 

[Eig96] M. Eigcn. Prionics or the kinetic basis of prion diseases. Biophys. Chem., 
63:A1-A18, 1996. 

[GPMW06] M. L. Greer, L. Pujo-Menjouet, and G. F. Webb. A mathematical analysis of 
the dynamics of prion proliferation. J. Theor. Biol., 242(3):598-606, Oct 2006. 

[Lau02] P. Laurcngot. The discrete coagulation equations with multiple fragmentation. 
Proc. Edinb. Math. Soc. (2), 45(l):67-82, 2002. 

[Len09] N. Lenuzza. In preparation. PhD thesis, Paris, 2009. 

[LM07] P. Laurcngot and S. Mischler. From the discrete to the continuous coagulation- 
fragmentation equations. Proceedings of the Royal Society of Edinburgh: Section 
A Mathematics, 132(05):1219-1248, 2007. 



23 



[LW07] 



P. Laurencot and C. Walker. Well-poscdness for a model of prion proliferation 
dynamics. Journal of Evolution Equations, 7(2):241-264, May 2007. 



[Mal82] P. Malliavin. Integration et probabilities. Analyse de Fourier et analyse spectrale. 
Masson, Paris, 1982. Collection: Maitrise de Mathcmatiques Pures. 

[MGA05] J. Masel, N. Genoud, and A. Aguzzi. Efficient inhibition of prion replication by 
PrPFc2 suggests that the prion is a PrPSc oligomer. J. Mol. Biol, 345:1243- 
1251, 2005. 

[MJN99] J. Masel, V. A. Jansen, and M. A. Nowak. Quantifying the kinetic parameters 
of prion replication. Biophys. Chem., 77(2-3):139-152, Mar 1999. 

[Per07] B. Pcrthame. Transport equations in biology. Frontiers in Mathematics. 
Birkhauser Verlag, Basel, 2007. 

[Rez08] H. Rezaei. Prion protein oligomerization. Current in Alzheimer research, 5:572- 
578, 2008. 

[SRH+05] J. R. Silvcira, G. J. Raymond, A. G. Hughson, R. E. Race, L. Sim, V, S. F. 

Hayes, and B. Caughcy. The most infectious prion protein particles. Nature, 
437:257-261, 2005. 



24 



